library(findMFDIA)

con=createConnection("/home/witold/AQUAYest/napedro_L120228_001_SW_0.sqlite")
addIndexMZRT2Features(con)
addIndexRT2Features(con)

sumAr=summaryDB(con)

createBestFeaturesView(con,rtext=0,mzext=0,Volume=0)
map16 <- getFeaturesFlex(con, 20)
class(map16)
class(con)
plot(map16$MZ,map16$RT,pch="+")

map16$rtProjection[[1]]
plot(readBin(map16$rtProjection[[1]],numeric(),n=100,size=4))

map16 <- getFeaturesFlex(con, 16, massrange=c(1000,1010),rtrange=c(400,8000))

